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P_i' Abstract 

In D-term inflation models, Affleck-Dine baryogenesis produces isocurvature 
density fluctuations. These can be perturbations in the baryon number, or, in 
the case where the present neutralino density comes directly from B-ball decay, 
perturbations in the number of dark matter neutralinos. The latter case results 
in a large enhancement of the isocurvature perturbation. The requirement that 
the deviation of the adiabatic perturbations from scale invariance due to the 
Affleck-Dine field is not too large then imposes a lower bound on the magnitude 
of the isocurvature fluctuation of about ICT 2 times the adiabatic perturbation. 
This should be observable by MAP and PLANCK. 
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The quantum fluctuations of the inflaton field give rise to fluctuations of the energy 
density which are adiabatic [|TJ. However, in the minimal supersymmetric standard 
model (MSSM), or its extensions, the inflaton is not the only fluctuating field. It is 
well known that the MSSM scalar field potential has many flat directions 0, along 
which a non-zero expectation value can form during inflation, leading to a condensate 
after inflation, the so-called Affleck-Dine (AD) condensate jj. The AD field is a com- 
plex field and, in the currently favoured D-term inflation models [|J on which we focus 
in this letter, is effectively massless during inflation. Therefore both its modulus and 
phase are subject to fluctuations. In D-term inflation models the phase of the AD field 
receives no order H corrections after inflation and so its fluctuations are unsuppressed 
|J. Because the subsequent evolution of the phase of the AD condensate generates 
the baryon asymmetry ||, the fluctuations of the phase correspond to fluctuations in 
the local baryon number density, or isocurvature fluctuations, while the fluctuations 
of the modulus give rise to adiabatic density fluctuations. We will show that the adia- 
batic fluctuations may in fact dominate over the inflaton fluctuations, with potentially 
adverse consequences for the scale invariance of the perturbation spectrum, thus im- 
posing an upper bound on the amplitude of the AD field. As a consequence, there is a 
lower bound on the isocurvature fluctuation amplitude. The magnitude of this lower 
bound will depend on the nature of the AD field. D-term inflation models require 
that d > 4, where d is the dimension of the non-renormalizable superpotential terms 
stabilizing the potential, in order to avoid thermalizing the AD field too early ||], 
whilst R-parity conservation, required to eliminate dangerous renormalizable B and L 
violating terms from the MSSM, rules out odd values of d. Therefore we will focus on 
the d = 6 direction, in particular the u c d c d c direction, in the following. We will show 
that the isocurvature fluctuations will most likely be observable in forthcoming cosmic 
microwave background (CMB) satellite experiments. 

An important point is that the AD condensate is not stable but typically breaks 
up into non-topological solitons || 0] which carry baryon (and/or lepton) number 
[§, [| and are therefore called B-balls (L-balls). This is a generic feature which is not 
realized only if the fluctuations take the AD field into certain leptonic d = 4 (" H U L" ) 
directions. The formation of the B-balls takes place with an efficiency likely to 
be in the range 0.1 to 1 fl0| . Hence the AD isocurvature fluctuations are inhereted 
by the B-balls. The properties of the B-balls depend on SUSY breaking and on the 
flat direction along which the AD condensate forms. We will consider SUSY breaking 
mediated to the observable sector by gravity. In this case the B-balls are unstable but 
long-lived, decaying well after the electroweak phase transition has taken place J7J, with 
a natural order of magnitude for decay temperature ~ 0(1) GeV. This assumes a 
reheating temperature after inflation, Tr, of the order of 1 GeV. Such a low value of 
Tr is necessary in D-term inflation models because the natural magnitude of the phase 



of the AD field, 5qp, is of the order of 1 in D-term inflation and along the d=6 direction 
AD baryogenesis implies that the baryon to entropy ratio is r\ B ~ 5cp(T R /10 9 GeV) 
[ P"T|j . It is significant that a low reheating temperature can naturally be achieved in 
D-term inflation models, as these have discrete symmetries in order to ensure the 
flatness of the inflaton potential which can simultaneuously lead to a suppression of 
the reheating temperature |TT] . 

Because the B-ball is essentially a squark condensate, in R-parity conserving mod- 
els its decay produces both baryons and neutralinos (x), which we assume to be the 
lightest supersymmetric particles (LSPs), with n x ~ 3n# [H], 0. This case is par- 
ticularly interesting, as the simultaneous production of baryons and neutralinos may 
help to explain the remarkable similarity of the baryon and dark matter neutralino 
number densities JTU], [Tj|. With B-ball decay temperatures Td ~ Oil) GeV, the decay 
products no longer thermalize completely and, so long as T<j is low enough that they do 
not annihilate after B-ball decay ||12|| , retain the form of the original AD isocurvature 
fluctuation. Therefore in this scenario the cold dark matter particles can have both 
isocurvature and adiabatic density fluctuations, resulting in an enhancement of the 
isocurvature contribution relative to the baryonic case. On the other hand, if the neu- 
tralinos from B-ball decays annihilate, the neutralino contribution to the isocurvature 
fluctuation will be erased, leaving only the baryonic contribution. Although we will 
be primarily interested in the neutralino isocurvature fluctuation case, we will also 
comment on the purely baryonic case in the following. 

Isocurvature perturbations have been studied previously in particular in the 
context of axion models JEJJ [15]. The isocurvature perturbations give rise to extra 
power at large angular scales but are damped at small angular scales. The amplitude 
of the rms mass fluctuations in an 8h~ l Mpc -1 sphere, denoted as eg, is about an order 
of magnitude lower than in the adiabatic case. Hence COBE normalization alone is 
sufficient to set a tight limit on the relative strength of the isocurvature amplitude. 
Small isocurvature fluctuations are, however, beneficial, in that they improve the fit 
to the power spectrum in Oq = 1 CDM models with a cosmological constant [[14] (or 
Qo = 1, A = CDM models with some hot dark matter For instance, in the 

context of axion models it has been found |14[ that an Q = 1 mixed fluctuation model 
with a relative isocurvature perturbation amplitude of 5%, Q a = 0.4 and = 0.6 
would give a very good fit to the data. However, the isocurvature fluctuations seem 
to require a large axion decay constant, which is already excluded unless there is 
considerable late entropy production fll4 |. The Affleck-Dine case we consider here is 
more economic, in the sense that it requires only the particles of the MSSM. 

In D-term inflation models, the AD field <3? = (f>e %e / \/2 = (<pi + z0 2 )/\/2 remains 
effectively massless during inflation. Therefore the real fields (pi are subject to quantum 



fluctuations with 

^(x) = Wj ^e"^ x 5 k , (1) 

where V is a normalizing volume and where the power spectrum is the same as for the 
inflaton field, 

where Hi is the value of the Hubble parameter during inflation. Thus, for given 
background values 9 and 0, (with 9 naturally of the order of 1) one finds 

( S9 \ Hj Hjk-W 

\Tan{9)J k Tan{9)<P y/2Tan{9)(j) I ' 

where <pi is the value of <fi when the perturbation leaves the horizon. After inflation, 
during the inflaton oscillation dominated period, the mass squared of the magnitude 
of the AD field will receive an order H 2 correction, which must be negative in order 
to have a non-zero </> and so AD baryogenesis 0, whilst its phase receives no order 
H corrections. Therefore, the magnitude of the AD field $ remains at the non-zero 
minimum of its potential until H ~ 777.5, where mg ~ 100 GeV is the SUSY breaking 
scalar mass, whence it begins to oscillate and the baryon asymmetry tx Sin(9) 
forms. Since 9 and 59 remain constant until H ~ rug, we have 
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with 59 /Tan{9) given by Eq. (|3|). 

We first consider the case where the adiabatic perturbation is mostly due to the 
inflaton. The adiabatic perturbation is determined by the invariant ( = 5p/(p + p) 
with 5p = V'5<f). During inflation, when all the fields are slow rolling, one finds |15| 



where <5 7 = 5p 7 /p 7 . For super-horizon size isocurvature fluctuations 5p/p = 0, so that 
m x 5n x + rriB5n B + 4(p 7 + p v )5T /T = (here p 7 and p v ~ 0.68p 7 are respectively 
the photon and the neutrino densities, and we assume for simplicity that there are no 
massive neutrinos). We then find that in the presence of both adiabatic (5(n x /s) = 0) 
and isocurvature (5(n x /s) 7^ 0) fluctuations 

5_T = p x 6<p + PB 5 ( S 

T 3 (p x + p B )+4(p 7 + Pi ,) ' 1 } 



where 5 X = 5n x /n x for non-relativistic particles x. The isocurvature fluctuations of the 



LSPs are related to the baryonic isocurvature fluctuations by 5n^ = 3/b^"-b\ with 



Sn^ given by Eq. (|J). In the linear perturbation theory adiabatic and isocurvature 
fluctuations evolve independently so that the total perturbation is just the sum of the 
two. 

The total LSP number density is the sum of the thermal relic density and the 
density nj?' = 3fBns originating from the B-ball decay. Using Eq. (|6]), the isocurva- 
ture fluctuation imposed on the CMB photons is then found to be 



7 




4 (1 + sgfe) P^S. 
3(P X + Pb) + 4(p 7 + Pu) 



B 




where p^ is the LSP mass density from the B-ball, Vt m (O x ) is total matter (LSP) 
density (in units of the critical density), and is given by Eq. @. To obtain the last 
line in Eq. (0), we have used the fact that p 7 is negligible. In the notation of reference 



15] and using Eq. (Bj) we can write 




where S is the inflaton field with a potential V(S) and M = M P i/y/87r. 

In the simplest D-term inflation model, the inflaton is coupled to the matter fields 
and ip + carrying opposite Fayet-Iliopoulos charges through a superpotential term 

W = nSip_'ip + 0, |[. At one loop level the inflaton potential reads 

where £ is the Fayet-Iliopoulos term and g the gauge coupling associated with it. 
COBE normalization fixes £ = 6.6 x 10 15 GeV ||16| . In addition, we must consider 
the contribution of the AD field to the adiabatic perturbation. During inflation, the 
potential of the d = 6 flat AD field is simply given by 

V{6) = -^-0 10 . (10) 

With p = V(S) + V((f>) and p + p = S 2 + 2 one finds, taking both S and to be slow 
rolling fields with S = -V'(S)/(3H I ) and = -K'(0)/( 3 #/)> that the invariant ( is 
now 

v\4>) + v\S) 

where we have used the fact that both fields are massless, so that 5S = 5(f). Thus 
the field which dominates the spectral index of the perturbation will be that with the 



largest value of V and V" . The index of the power spectrum is given by n = l+2r/ — 6e, 
where e and rj are defined as 

, = . (12, 

The present lower bounds imply |An| < 0.2. (This bound will be much improved by 
future satellite missions). In the case where the derivatives with respect to the inflaton 
dominate (for which the potential is dominated by Vq for all £ < M), \An\ = l/N « 
0.02 for ~ 50. Once the derivatives with respect to the AD field dominate, the 
spectral index increases rapidly with <p\ from rj (e), | An\ is proportional to <j) 8 (</> 18 ). The 
condition for the AD field to dominate the spectral index is that <fi > Max(<p Cl , C2 ), 
where 

0.64(# 3 A~ 2 £ 4 M 5 ) 1/9 (13) 



0.48 ( f V /4 (M£) 1/2 , (14) 



and 

rh. ~ D.4S ( 

,A 

and where we have used the fact that during the slow roll-over S"^ ~ g 2 M 2 N/(47r 2 ). 
The inflaton derivatives will dominate once <fi < Min((j) Cll (j) C2 ). (In practice C1 and 
4> C2 only differ by a factor of less than 2). As a result of the rapid increase of the 
spectral index once the AD derivatives dominate, the condition that the spectral index 
is acceptably close to scale invariance essentially reduces to the condition that it is 
dominated by the inflaton. The lower bounds on (3 corresponding to C1 and C2 are 
then 

P > C1 ~ 6.5 x 10-V /3 A 4/9 ^ 2 Tan(^)- 2 (15) 

and 

P > (3 C2 = 2.5 x 10-V /2 A 1/2 ^ 2 Tan(^)- 2 . (16) 

(For most values of the couplings the latter leads to a slightly more stringent lower 
bound). Thus significant isocurvature fluctuations are a definite prediction of the AD 
mechanism. 

There are two limiting cases: /# ^> m B /?>f B m x , for which u = 1, and fs ^ 
TnB/3fBm x > f° r which ^ = ^s/^m- The latter corresponds to the case where the 
B-balls form very inefficiently or where the neutralino contribution to the isocurva- 



ture perturbation is erased by annihilations [It], |12|1 . In this case only the baryonic 
isocurvature perturbation remains. 

The actual lower limit on P depends on the unknown couplings g and A, as well 
as on 9. To obtain an estimate for /3 C2 in the case where the dark matter neutralinos 
come directly from B-ball decay, let us adopt the following values: g ~ ^gut — 
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Figure 1: The relative difference ACj/Cj between the purely adiabatic and a mixture 
of adiabatic and isocurvature angular power spectra with (3 = 0.001 (dotted line) and 
(3 = 0.0001 (solid line) for a purely CDM Q = l model (with Q B = 0.05, h = 0.5 and 
the spectral index n — 1). Shown is also the projected PLANCK error level, averaged 
over ten multipoles (dashed line). 



0.7, A ~ 1/5! = 0.008 (corresponding to non-renormalizible interaction with physical 
strength set by M || PH), Tan{9) 2 ~ 1 (corresponding to 9 « 7r/4) and, assuming 
/b is not too small, u; ~ 1. We then find that (3 > (3 C2 ~ 1.3 x 10 -3 , with the 
conservative lower bound perhaps an order of magnitude smaller, (3 ~ 10~ 4 . For the 
purely baryonic case, the value of u depends on fl B and fl m . Nucleosynthesis combined 
with the current best estimate of the expansion rate (0.6 < h < 0.87 [17]]) implies that 

0. 006 < Qb ^ 0.036. Thus with fi m = 1 (0.4) we obtain that u for neutralino case is 
30-150 (10-60) times larger than in the purely baryonic case. Therefore in the baryonic 
case the corresponding value of j3 is two to four orders of magnitude smaller. Thus 
there is a significant enhancement of the isocurvature perturbation in the case where 
the dark matter neutralinos come from B-ball decay without subsequent annihilations. 
It should be emphasized that there is no physical reason to expect to be close to its 
upper bound, so (3 may, in general, be expected to be much larger that these lower 
bounds. In this case, even without the neutralino enhancement, the purely baryonic 
isocuvature fluctuation may well be important. 

In Fig. 1 we display the difference between the purely adiabatic power spectrum 
and the spectra with (3 ^ 0. We also plot the expected error for the Planck Surveyor 
Mission, following the estimates in ref. |18|]. (A similar error is expected for MAP 
for I < 500). The standard error reads (AC,) 2 = 2(Q + 5) 2 /[(2l + l)/ sky ], where / sky 
is the fraction of the sky sampled (we take / s k y = 0.65) and S is from the beam, the 
angular resolution and the sensitivity, as discussed in jl8fl ; 5 becomes non-negligible 
only for / > 1000 for PLANCK and / > 500 for MAP. One should bear in mind that, 
in principle, each multipole provides an independent measurement of the spectrum. 
As can be seen, detecting isocurvature fluctuations at the level of (3 ~ 10~ 4 should 
be quite realistic by averaging over a sufficient number of multipole measurements. 
However, detecting or setting an actual lower limit on f3 will require a much more 
careful analysis, which we do not attempt here. Nevertheless, on the basis of Fig. 

1, it seems likely that the forthcoming CMB experiments will definitely be able to 
see isocurvature perturbations in the case where the baryons and neutralinos come 
directly from the decay of unstable B-balls in the context of D-term inflation models, 
hence offering a test not only of the inflationary Universe but also of the B-ball variant 
of AD baryogenesis. 

In conclusion, AD baryogenesis in the context of D-term inflation generally implies 
the existence of isocurvature density fluctuations. In the case where B-balls, which 
are generally expected to form in AD baryogenesis, decay late enough to produce the 
observed neutralinos without annihilations, the isocurvature fluctuctions should be 
observable by MAP and PLANCK. Even in the case where only baryonic isocurvature 
fluctuations arise, there is still a reasonable possibility of observing them, although 
in this case it is less certain. Thus isocurvature fluctuations are a clear fingerprint 



of D-term inflation and Affleck-Dine baryogenesis. In particular, for the case where 
the neutralino dark matter comes directly from B-ball decays, which allows for an 
understanding of the remarkable similarity of the baryon and dark matter number 



densities [pT0| , fL2|| , observation of isocurvature perturbations combined with the non- 
thermal nature of the dark matter neutralino density (testable by observation of the 
sparticle spectrum |12|]) would strongly support the late decaying B-ball scenario and 
D-term inflation, giving us a deep insight into the nature of particle physics and the 
very early Universe. 
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